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ABSTRACT 

In previous work, we identified six Sun-like stars observed by Kepler with exceptionally clear asteroseismic signatures of rotation. 
Here, we show that five of these stars exhibit surface variability suitable for measuring rotation. We compare the rotation periods 
obtained from light-curve variability with those from asteroseismology in order to further constrain differential rotation. The two 
rotation measurement methods are found to agree within uncertainties, suggesting that radial differential rotation is weak, as is 
the case for the Sun. Furthermore, we find significant discrepancies between ages from asteroseismology and from three different 
gyrochronology relations, implying that stellar age estimation is problematic even for Sun-like stars. 
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1. Introduction 


Recently, asteroseismology has become a valuable tool to study 
the internal rotation of stars. This has been done on a variety of 


different stars (see, e.g., |Aerts et al.|20 03; Charpinet et al. 


Kurtz et al. 


2014|), including Sun-like stars ( [Gizon et al. 


2009 


2013 


Davi es et al. 12015] ). Asteroseismic measurements of rotation in 


Sun-like stars have generally been limited to the average inter¬ 
nal rotation period. More recently, Nielsen et al.] |2014| ) identi¬ 
fied six Sun-like stars where it was possible to measure rotation 
for independent sets of oscillations modes and found that radial 
differential rotation is likely to be small. 

From the Sun, we know that different methods for measur¬ 
ing rotation, e.g., spot tracing (|D’ Silva & Howard) 1994| and he¬ 
lioseismology ( |Schou et al.| 1998| ), show the outer envelope of 
the Sun rotating differentially. The solar surface rotation period 
changes by approximately 2.2 days between the equator and 40° 
latitude ( [Snodgrass & Ulrich| 199Q| ), while in the radial direction 
the rotation period changes by ~ 1 day in the outer few percent 
of the solar radius ( |Beck|2000| . While small on the Sun, differ¬ 
ential rotation could be larger on othe r stars . With data available 
from the Kepler mission ( [Borucki et al.|2010| , it is now possible 
to compare the measurements of rotation using both asteroseis¬ 
mology and observations of surface features. 

Differences between the results of the two methods would 
imply that differential rotation exists between the regions where 
the two methods are sensitive; namely, the near-surface layers for 
asteroseismology, and the currently unknown anchoring depth 
of active regions in Sun-like stars. However, given the similar 
scales of the radial and latitudinal differential rotation seen in 
the Sun, both of these mechanisms may contribute to an ob¬ 
served difference. On the other hand, agreement between the 
two methods would suggest weak differential rotation. Further¬ 
more, this would also mean that we can measure rotation with 
asteroseismology in inactive stars that exhibit few or no surface 
features at all, thus giving us an additional means of calibrating 


gyroc hronology relations (e.g., Skumanich 1972} |Barnes||2007 


2010 ), which so far have primarily relied on measuring the sur¬ 
face variability of active stars. 

However, the acoustic modes (or /7-modes) used to measure 
rotation in Sun-like stars are damped by surface magnetic ac¬ 
tivity ( [Chaplin et al.||2011| ). It is therefore difficult to find stars 
with both a high signal-to-noise oscillation spectrum and coher¬ 


ent signatures of rotation from surface features. Nielsen et al. 


( [2014 ) identified six Sun-like stars from the Kepler catalog with 
exceptional signal to noise around the p-mode envelope. Five 
of these stars exhibit rotational variability from surface features. 
We compare the surface variability periods with those measured 
from asteroseismology, with the aim of further constraining the 
near-surface differential rotation of these stars. 


2. Measuring rotation periods 

We use the ~58 second cadence observations from the Kepler 
satellite to capture the high frequency acoustic oscillations, and 
the ~29 minute cadence observations to search for surface vari¬ 
ability. The stars we identified were chosen because of their sim¬ 
ilarity to the Sun with respect to surface gravity and effective 
temperature, and for having a high signal-to-noise ratio at the p- 
mode envelope. We use all the available dat^for each star, from 
quarter 1 to 17, spanning approximately four years. 

2.1. Asteroseismic rotation periods 

The oscillation modes in a star can be described with spherical 
harmonic functions, characterized by the angular degree / and 
azimuthal order m, and a radial component with order n. When 
the star rotates the modes with azimuthal order m > 0 become 

1 Downloaded from the Mikulski Archive for Space Telescopes, 
http: //archive, stsci. edu/kepler/ 
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Fig. 1. Power density spectrum of the short cadence white light pho¬ 
tometry time series of KIC005184732. Gray represents the spectrum 
smoothed with a 0.1/zHz wide Gaussian kernel, with the best-fit model 
shown in red. Black dashed indicates the surface variability frequency. 
The inset shows a zoom of the p-mode envelope at an l = 2,0 pair, 
where the black markers indicate the individual m-components. The 
separation between components of the l - 2 mode is the rotational split¬ 
ting. 


Doppler shifted, and if this can be measured, the rotation of the 
star can be inferred. For Sun-like stars, the oscillation modes 
have a sensitivity weighted toward the surface of the star, and so 
the asteroseismic measurements predominantly probe the rota¬ 
tion of the star in this region (see, e.g.,[Lund et al. |2014| ). 

To measure this frequency splitting we fit a model to the 
power spectral density. The modes in a solar-like oscillator are 
stochastically excited and damped oscillations can be described 
by a series of Lorentzian profiles, 




&\m (0 $ nl 


1=0 m=-l 


+ (2/T n i m ) 2 (v Vnlm) 


+ *(v), 


( 1 ) 


where v is the frequency, and B(v) is the background noise. Here 
the sums are over the radial orders (typically -8 values of n ), 
the angular degrees l < 3 and the azimuthal orders -l < m < l. 
Not all radial orders show clear / = 3 modes, and so only a 
few l = 3 modes are included in the fits, depending on the star. 
The mode power in a multiplet, S n u is a free parameter. The 
mode visibility $i m (i) is a function of the inclination angle i of 
the rotation axis relative to the line of sight (see Gizon & Solanki 
2003). We assume here that all modes share the same value of i. 

The frequencies v„/ m of the modes reveal the rotation infor¬ 
mation. For Sun-like stars rotating no more than a few times the 
solar rotation rate the mode frequencies can be parametrized as 
Vnim - v„/+mbv, where Sv is the rotational splitting. In our model 
fit we assume that the stellar interior rotates as a solid body, so 
that all modes share a common rotational splitting. 

The full widths Y n i m = T(v n i m ) of the Lorentzian profiles de¬ 
pend on the mode lifetimes, which in turn depend on the mode 
frequencies v„/ m . In the Sun this variation with frequency can be 
modeled by a low-order polynomial ( |Stahn|20l0| ). We therefore 
opt to parametrize the mode widths as a third-order polynomial 
in (v - v max ), where v max is the frequency at maximum power of 
the p-mode envelope. Note that only the / = 0 modes contribute 
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Fig. 2. Peak periods in the Lomb-Scargle periodogram as a function of 
observation quarter for KIC004914923 (top) and KIC005184732 (bot¬ 
tom). Solid gray shows the median period, and dashed gray shows the 
median absolute deviation. 


to the fit of this function since they are not broadened or split by 
rotation. 

The background noise level B (v) is caused by several terms: 
the very long-term variability (hours to days) from magnetic 
activity, the variability (tens of minutes) from stellar surface 
granulation, and the photon noise. For asteroseismology appli¬ 
cations, these noise components can be adequately described by 
two Harvey-like background terms (e.g., |Handberg & Campante] 
2011) and a constant term to account for the white photon noise. 


We find the best fit with maximum likelihood estimation us 
ing a Markov Chain Monte Carlo sampler (Foreman-Mackey 


|et al.|2013] ). The best-fit values of each parameter are estimated 
by the median of the corresponding marginalized posterior dis¬ 
tribution, and the 16 th and 84 th percentile values of the distri¬ 
butions represent the errors. We fit the entire spectrum, includ¬ 
ing all mode parameters and background terms simultaneously 
(a so-called global fit), in contrast to [Nielsen et al.| ( [2014| ) who 
fit the background and separate sets of p-modes individually. A 
comparison between the global rotational splitting obtained in 
this work, and the mean of the rotational splittings measured by 
|Nielsen et al.| ( |2014) are shown in Table [T] 

An example spectrum of KIC005184732 is shown in Fig. |T] 
The background terms dominate the low frequency end of the 
spectrum, while the p-mode envelope appears clearly above the 
noise level at -1800 p Hz. The inset shows a part of the p-mode 
envelope, where the splitting of an l = 2 mode is clearly visible. 


2.2. Surface variability periods 

The rotation of stars can also be inferred from the periodic vari¬ 
ability of their light curves caused by surface features on the 
stellar disk, such as active regions. For the Sun, active regions 
are good tracers of the surface rotation of the plasma, to within 
a few percent ( |Beck||2000| ). In other stars this difference could 
potentially be larger. 

From the periodic variation of the stellar light curve it is rel¬ 
atively straightforward to identify the rotation period of the star, 
and this has been done using automated routines for tens of thou¬ 
sands of stars ((Nielsen et"aL]|2013t [Reinhold & Reiners|2013t 
|McQuillan et al.||2014| ). However, these routines typically rely 
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on coherent surface variability over long time periods. The stars 
studied here do not show simple and regular surface variability 
and, therefore, they do not appear in these catalogs. Furthermore, 
the PDC_MAP pipeline ( [Smith et al.|2012[|Stumpe et al.|2014 ) 
is known to suppress variability on timescales longer than ~20 
days ^Christiansen et al .|]2013| ) 

We therefore manually reduce the raw pixel data for each 
star, and search this and the PDC_MAP reduced data for signs 
of rotational variability. To reduce the raw data, we use the kep- 
cotrend procedure in the PyKe software packag<0 The PyKe 
software uses a series of so-called cotrending basis vectors 
(CBVs) in an attempt to remove instrumental variability from 
the raw light curves. The CBVs are computed based on vari¬ 
ability common to a large sample of stars on the detector, and 
so in principle represent the systematic variability. The number 
of CBVs to use in the reduction is not clear and we therefore 
compute several sets of reduced light curves using a sequentially 
increasing number of CBVs (up to six) for each star. We require 
that the variability appears in light curves reduced with different 
numbers of CBVs. This minimizes the risk of variability signa¬ 
tures that are caused by the reduction. 

To measure the rotation period of each star we used the 
method of Nielsen ~et al.|(|2013|). We compute the Lomb-Scargle 
periodogram ( |Lomb[ 1976[ |Sca rgle 1982 ) of each quarter and 
identify the peak of maximum power for periods less than 45 
days. The median period of these peaks is used as a first order es¬ 
timate of the rotation period. We then assume that active regions 
appear at similar latitudes, and so should have similar periods. 
This is done by requiring that peaks must lie within 4 median 
absolute deviations (MAD) from the median. The median and 
MAD of the remaining peak periods are the measured rotation 
period and error on the rotation period. Examples are shown in 
Fig. [2] The rotation periods of th ese stars measured using this 
method agree with those found by Garci a et aT7| ( |2Q 1 4| ) |do Nasci- 
|mento et~aL] ( |2Q14| ). 


3. Comparing asteroseismic and surface variability 
periods 

In Table [T] we list the periods for each analysis. These values 
are also shown in Fig. [3] The asteroseismic and surface vari¬ 
ability periods agree very well for all the stars analyzed here. 
However, the asteroseismic rotation period measurements are 
systematically lower than those from surface variability. As a 
first order significance estimate of this difference we fit a lin¬ 
ear function P s = aP a to the measured values. We find the best 
fit by maximum likelihood estimation, using the marginalized 
posterior distributions for the asteroseismic measurements, and 
assume Gaussian posterior distributions for the surface variabil¬ 
ity periods. The best-fit solution returns a = 1.044^ ^, show¬ 
ing that the difference between the two measurement methods 
is insignificant (see Fig. [3]). As a consistency check we perform 
an identical fit to the average rotation periods P n from |Nielsen| 
|et al.| ( |2014| ), assuming Gaussian errors (see Table [T}. We find 
a slope of 1.006 ^-^q, in good agreement with the slope of Pa- 
However, this does not account for the non-Gaussian form of the 
marginalized posteriors of the fits performed in that work. 


4. Gyrochronology 

Here we estimate the stellar ages from three different gy¬ 
rochronology relations provided by [Mamajek & Hillenbrand) 

2 http://keplerscience.arc.nasa.gov/PyKE.shtml 



Fig. 3. Surface variability period P s as a function of asteroseismic pe¬ 
riod Pa- The 1:1 line is shown with a black dashed line. The best-fit 
value and errors to the slope are shown with red solid and dashed lines, 
respectively. For comparison the Carrington period of 27.3 days is also 
shown (this is not included in the fit). 
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Fig. 4. Differences between gyrochronology ages and asterosei smic 
ages for three different gyrochronology relations. The relations are|Ma-| 
|majek & Hillenbrand| ( |2008| ), |Barnes < |2010| ), and |Garcfa et al.| ( [2014| ). 


( 2008]), [Barnes| ( |2010| ), and Gar cia et aL] ( |2014| ). These relations 
represent distinct methods of estimating gyrochronology ages. 
In all these relations, we use the asteroseismic rotation periods 
from Sect. 12.11 


The type of gyro-relation used by Mamajek & Hillenbrand 
( 2008[ ), first proposed byjBarnes| (2007| , assumes that the rota¬ 
tion period of a star varies as powers of age and mass, where 
the B - V color is used as a proxy for the latter. This relation is 
calibrated to young cluster stars with ages < 625 Gyr, with only 
the Sun representing the older stellar populations. For the stars 
studied here we use B - V values from |H0g et al.|(|2000 ), using 
extinction values from the KIC catalog (Brown et a l.|2011 ). 

The relation by Garci a et al.| ( |2014| ) is calibrated to stars with 
asteroseismically determined ages from |Mathur et al.| ( |2012[ ) and 
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Table 1. List of measured and computed parameters for the five stars. 


KIC004914923 KIC005184732 KIC006116048 KIC006933899 KIC010963065 

Mass [M©] 

Radius [R©] 

X c 

r [days] 

T eff [K] 

P s [days] 

Pa [days] 

P n [days] 

Seismic age [Gyr] 


Gyrochronolo gy ages [Gyr] 


Barnes 

(201C 

l) 




2-68!^ 

9 79+O.47 
' -0.54 

3 09 +0 - 36 
-0.44 

5.58!^ 

1 Q7+0.49 

L y ' -0.34 

Mamajek & 1 

Tillenbrand 

(2008 

) 

3.56!j^ 

2 36 +0,63 

^-0.53 

3 36 +1 14 
J,JU -0.82 

O C7+3.73 
-3.08 

4 38 +3 - 77 

H.JO_i 74 

Garcia et al. ( 

2014 

) 



3-25!^ 

3 65 +1 - 37 

-0.95 

3 05 +0,92 

-0.67 

7 09+4.78 
/. 53-2 84 

1 4Q+0.51 
-0.34 


1.118 + 0.020 

1.205 + 0.025 

1.023 + 0.021 

1.096 + 0.026 

1.062 + 0.021 

1.378 + 0.009 

1.342 + 0.010 

1.225 + 0.008 

1.574 + 0.025 

1.220 + 0.009 

0.000 

0.216 

0.036 

0.000 

0.108 

16.94 

23.30 

34.34 

28.25 

33.03 

5880 + 70 

5865 + 70 

5990 + 70 

5870 + 70 

6090 + 70 

19.49 + 3.12 

20.69 + 0.50 

17.96 + 2.11 

31.63 + 1.43 

12.27 + 0.32 

17.98! 3 ;* 7 

IQ 44+ 1 - 63 

17 61 +0 - 95 

1 /,u -1.31 

29.92+ 4 - 9 ° 

12 01 +L42 

22.17 + 3.14 

18.00+1.76 

16.34 + 0.89 

28.65 + 5.53 

14.45 + 1.43 

6.23 + 0.36 

4.39 + 0.13 

5.70 + 0.21 

6.57 + 0.30 

4.18 + 0.19 


Notes. The stellar masses and radii, core hydrogen fraction X c , and convective turnover time r, as well as the seismic ages are computed from 
the best-fit models as described in Section [4] The effective temperatures are adapted from |Bruntt et al.| ( [2012| ). The measured surface variability 
rotation periods P s and asteroseismic rotation periods P A are also listed. These periods are compared to the averaged asteroseismic rotation periods 
P N from|Nielsen et al.|{|2014|). The P N values are a verages over several sets of oscillation modes. The gyrochronology ages are calculated from 
|Bames|p010| ), |Mamajek & Hillenbrand| ( |2008| ), and |Garcfa et al.| ( [20T4] >. 


[Metcalfe et al.||MT4| ). However, this relation does not include 
a mass dependence since the calibrators are all similar in mass 
(-1M©), and so this relation only depends on the rotation period 
of the star. 


Lastly, the relation by Barnesj p010| incorporates the effects 
of a rapid spin-down phase, experienced by very young stars 
(-100 Myr), into the prediction of the stellar age. This relation 
uses convective turnover time r as a proxy for mass, and is also 
calibrated using young cluster stars, with the Sun as the only 
anchoring point for older stars. 

We determine the convective turnover times from stellar 
models fit to the stellar oscillation frequencies (see online ma¬ 
terial in [Nielsen et al.| ( |2014| )). These stellar models are pro¬ 
duced using the Modules for Experiments in Stellar Astro¬ 
physics (MESA, Paxton et al. 201 1[ |2013| ), with the same in¬ 
put physics as in |B all & Gizon| (|2014|). Initia l fits were deter¬ 
mined using the SEEK method (Quirion et al.|2010| ) to compare 
a grid of models to observed spectroscopic and global asteroseis¬ 
mic parameters. The models cover the mass range 0.60-1.60M o 
and initial metal abundances Z in the range 0.001-0.040. The 
helium abundance was presumed to follow the enrichment law 
Y = 0.245 + 1.45Z and the mixing length parameter was fixed at 
a Q = 1.908. 

As initial guesses, we used the median values of the grid- 
based fit and also generated ten random, uniformly-distributed 
realizations of the parameters. The 11 initial guesses were op¬ 
timized in five parameters (age, mass, metallicity, helium abun¬ 
dance, and mixing length) using a downhil l simplex ([N elder & 
Mead] [1965] ) to match spectroscopic data (Brun tt et al.||2012) 
and individual oscillation fr equencies, computed with ADIPLS 
(Christensen-Dalsgaard 200 8]) and corrected according to the cu¬ 
bic correction by Ball & Gizon| ( |2014| ). As in |Ball & Gizon| 
( |2014| ), the seismic and nonseismic observations are weighted 
only by their respective measurement uncertainties. Best-fit pa¬ 
rameters and uncertainties are estimated from ellipses bounding 
surfaces of constant^ 2 for all the models determined during the 
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optimizations, corresponding to a total of about 3000 models for 
each star. 

The best-fit models yield local convective velocities and we 
use these to compute r as in |Kim &~D emarque|( |1996| ), these are 
listed in Table Q] From these values of r and the asteroseismic 
rotation periods from section [271] we compute the gyrochronol¬ 
ogy ages using Eq. 32 from [Barnes| ( |2010) and the suggested 
calibration values therein. 

In addition to providing values of r, the stellar models also 
produce estimates of the stellar ages. In Fig. [4] these ages are 
compared with those derived from the three gyrochronology re¬ 
lations. It is immediately obvious that the gyro-relations are in 
general not consistent with the seismic ages. For some stars a 
particular gyro-relation may agree well with the seismic age, 
but for other stars that same relation deviates significantly, as 
seen, e.g., for the ages fro m |Garcia et al.| ( |2014| ) and |Ma- 
|majek & Hillenbrand ( |2008| ) relations for KIC005184732 and 
KIC010963065. However, it appears that the gyro-relations tend 
to underestimate the stellar ages relative to the seismic values. 

For the star KIC006933963, the gyro-relations predict an age 
older than the stellar models, which is likely because this star 
is a subgiant, with a slightly contracted core and expanded en¬ 
velope. Thus, by angular momentum conservation one expects 
the surface and near-surface regions to have slowed down faster 
than the gyro-relations predict, leading to an overestimated age. 
KIC004914923 also appears to be a subgiant based on its core 
hydrogen content (see Table [T]), but close inspection of the best- 
fit stellar model shows that this star has only just moved off the 
main sequence, and so may still adhere to these gyro-relations. 

It should be noted that the errors on the asteroseismic ages 
are internal errors from the fit. However, even conservatively as¬ 
suming a precision of -20% on the asteroseismic ages ([Aerts 
[2015| , the two methods of age estimation are still in poor agree¬ 
ment. 
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5. Conclusions 


We measured the rotation periods of five Sun-like stars, using as- 
teroseismology and surface variability. The measurements show 
that the asteroseismic rotation periods are ~4% lower than those 
from surface variability, but this was not found to be a statisti¬ 
cally significant difference given the measurement uncertainties. 
Thus there is no evidence for any radial (or latitudinal) differen¬ 
tial rotation between the regions that the two methods are sensi¬ 
tive to. 

Previous efforts to calibrate gyrochronology relations have 
focused on measuring rotation from surface variability of clus¬ 
ter stars. Given that rotation periods can be reliably determined 
by asteroseismology compared to the surface variability, gy¬ 
rochronology may also be applied to very inactive stars for 
which surface variability rotation periods are not available. 

Our comparison of three different gyrochronology relations 
from the literature found that these ages differed by up to sev¬ 
eral Gyr from the seismic ages. In general, the gyro-relations 
seemed to underestimate the stellar ages relative to the astero¬ 
seismic ages. It also demonstrated that the gyro-relations are not 
internally consistent. 

When considering these discrepancies it is important to re¬ 
member the calibrations used for each method. The iBarnesI 
( |2010| ) and [Mamajek & Hillenbrand ( |2008| ) relations are both 
calibrated using ages and rotation periods of the Sun and cluster 
stars, and so may not accurately predict ages of stars similar to or 
older than the Sun (e.g., [Meibom et al.|2015| ). On the other hand, 
asteroseismic ages from stellar models are only calibrated to the 
Sun, and in addition are dependent on the input physics and as¬ 
sumptions of the stellar models ( [Lebreton & Goupil|2014| . The 
relation by Garci a et al.| ( |2014| ) lacks a mass dependent term, in 
contrast to the two other tested relations. However, this was not 
expected to have a remarkable effect since the stars in our sample 
were approximately 1M 0 stars, similar to the calibrators used by 
Garcia et alT|( |2Q14| ). Despite this, there is still a discrepancy of 
several Gyr between the ages derived from this relation and our 
seismic ages for some of the stars. 

Tighter asteroseismic constraints for stellar modeling will 
become available when calibrations are obtained from clus¬ 
ter stars to be observed by space missions like K2 ( |Chaplin| 
et al.|2015] ) and the Planetary Transits and Oscillations of stars 


(PLATO) mission ( [Rauer et al.[[2014| . This in turn will allow 
us to reliably use asteroseismic rotation periods to further refine 
gyrochronology relations for magnetically inactive field stars. 
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